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Abstract 

It is shown that measuring the lepton invariant mass spectrum in inclusive 
semileptonic B X^H-v decay yields a model independent determination of 
\Vub\- Unlike the lepton energy and hadronic invariant mass spectra, nonper- 
turbative effects are only important in the resonance region, and play a para- 
metrically suppressed role when dF/dg^ is integrated over > {niB — rnD)"^ , 
which is required to eliminate the charm background. Perturbative and non- 
perturbative corrections are presented to order a^Po and Aqf^j-j/m^ , and the 
^qcd/^I corrections are used to estimate the uncertainty in our results. The 
utility of the B — > Xsi~^£~ decay rate above the ip{2S) resonance is discussed. 
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A precise and model independent determination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element V^b is important for testing the Standard Model at B factories via the 
comparison of the angles and the sides of the unitarity triangle. The first extraction of \Vub\ 
from experimental data relied on a study of the lepton energy spectrum in inclusive charmless 
semileptonic B decay [|1|. Recently \Vub\ was also measured from exclusive semileptonic 
B piu and B Tfiu decay 0, and from inclusive decays using the reconstruction of the 
invariant mass of the hadronic final state 

These determinations suffer from large model dependence. The exclusive \Vub\ measure- 
ments rely on form factor models or quenched lattice calculations at the present time.Q 
Inclusive B decay rates are currently on a better theoretical footing, since they can be com- 
puted model independently in a series in Aqcd/""^?) and as{mb) using an operator product 
expansion (OPE) [[7|-p!0[]. However, the predictions of the OPE are only model independent 
for sufficiently inclusive observables, while the B — > X^^z/ decay rate can only be measured 
by imposing severe cuts on the phase space to eliminate the ~ 100 times larger B — >■ X^iP 
background. For both the charged lepton and hadronic invariant mass spectra, these cuts 
spoil the convergence of the OPE, and the most singular terms must be resummed into 



a nonperturbative b quark distribution function ||TT|. While it may be possible to extract 



this from the photon spectrum in B ^ Xg'y |TT|JT^, it would clearly be simpler to find an 



observable for which the OPE did not break down in the region of phase space free from 
charm background. In this Letter we show that this is the situation for the lepton invariant 
mass spectrum. 

At leading order in the AqcD/mb expansion the B meson decay rate is equal to the b 
quark decay rate. Nonperturbative effects are suppressed by at least two powers of Aqcd/''^6- 
Corrections of order AQQY^/ml are characterized by two heavy quark effective theory (HQET) 



matrix elements P|-pi][, which are defined by 

X, = {B{v)\h^:^ {iDfh^\B{v))/2mB. 

A2 = {B{v)\^h^^a^M^-h^'> \B{v))/QmB. (1) 

These matrix elements also occur in the expansion of the B and B* masses in powers of 
AqcD/mfe, 

-r Al + 3A2 V Al — A2 /_v 

niB = mb + A h . . . , rriB* = nib + A h . . . . (2) 

znib 2,mb 

Similar formulae hold for the D and D* masses. The parameters A and Ai are independent 
of the heavy b quark mass, while there is a weak logarithmic scale dependence in A2. The 
measured B* — B mass splitting fixes A2(mb) = 0.12 GeV^. At C(AQQj-,/m^) seven additional 



A model independent determination of \Vub\ from exclusive decays is possible without first order 
heavy quark symmetry or chiral symmetry breaking corrections Q, but it requires data on ^ — >■ 
K*i~^i~. A model independent extraction is also possible from decays to wrong-sign charm 
but this is very challenging experimentally. See also for a discussion of extracting \Vub\ from a 
comparison of photon spectra in B and D radiative leptonic decays. 
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parameters arise in the OPE |T3|-|T^], and varying these parameters is often used to estimate 
the theoretical uncertainty in the OPE |14- |T^] - 

In inclusive semileptonic B decay, for a particular hadronic final state X, the max- 



imum lepton energy is 



(max) 



[mi 



m 



X 



)/2mB (in the B rest frame), so to elim- 



inate charm background one must impose a cut Ei > 



[ml 



— m 



D 



)/2mB- The maxi- 



mum lepton energy in semileptonic b quark decay is mb/2, which is less than the phys- 
ical endpoint mB/2. Their difference, A/2, is comparable in size to the endpoint region 
^^^(ondpoint) _ ^2^^2mB ~ 0.33 GcV. The effects which extend the lepton spectrum be- 
yond its partonic endpoint appear as singular terms in the prediction for dT/dEi involving 
derivatives of delta functions, ^'^"'^(l — 2Ei/mi,). The lepton spectrum must be smeared over 
a region of energies A£^^ near the endpoint before theory can be compared with experiment. 
If the smearing region Ai?^ is much smaller than Aqcd, then higher dimension operators in 
the OPE become successively more important and the OPE is not useful for describing the 
lepton energy spectrum. For AE^ ^ Aqcd? higher dimension operators become successively 
less important and a useful prediction for the lepton spectrum can be made using the first 



few terms in the OPE. When ^Ee ~ A, 



QCD, 



there is an infinite series of terms in the OPE 



which are all equally important. Since ^^^('^'^'1p°™*) jg about Aqcd, it seems unlikely that 
predictions based on a few low dimension operators in the OPE can successfully determine 
the lepton spectrum in this region. 



It was shown in |11| that the leading singularities in the OPE may be resummed into 
a nonperturbative light-cone distribution function for the heavy quark. To leading 

order in 1/mb, the effects of the distribution function may be included by replacing mt by 



m. 



= mf, + and integrating over the light-cone momentum 



dP 

dE, 



dTr 



dEf, 



(3) 



where dTp/dEi is the parton-level spectrum. Analogous formulae hold for other differential 
distributions [0,0. For purposes of illustration, we will use a simple model for the structure 
function given by the one-parameter ansatz []TU| 



32 



TT- 



'■A 



(1 — x)^ exp 



TT 



0(1 -x) 



X = 



A 



(4) 



taking the model parameter A = 0.48 GeV, corresponding to = 4.8 GeV. 

The charm background can also be eliminated by reconstructing the invariant mass 
of the hadronic final state, mx, since decays with mx < must arise from b —>■ u 
transition. While this analysis is challenging experimentally, the mx < rrin cut allows a 
much larger fraction oib ^ u decays than the E^ > {m\ — m|))/2ms constraint. This is 
expected to result in a reduction of the theoretical uncertainties pO| , pl| , although both the 
lepton endpoint region, E(^ > {m\ — m|))/2mB, and the low hadronic invariant mass region, 
mx < fTiD, receive contributions from the same set of hadronic final states (but with very 
different weights). However, the same nonperturbative effects which lead to the breakdown 
of predictive power in the lepton endpoint region also give large uncertainties in the hadron 
mass spectrum over the range ~ Am;, |]2^. In other words, nonperturbative effects 
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FIG. 1. The lepton energy and hadronic invariant mass spectra. The dashed curves are the 
h quark decay results to 0(as), while the solid curves are obtained by smearing with the model 
distribution function /(A;+) in Eq. (Q). The unshaded side of the vertical lines indicate the region 
free from charm background. The area under each curve has been normalized to one. 



yield formally 0{\) uncertainties in both cases, because numerically m|, ~ Aqcd^b- The 
situation is illustrated in Fig. |l|. 

The situation is very different for the lepton invariant mass spectrum. Decays with 

= [pi + > {iTT-B — iTioY must arise from b u transition. Such a cut forbids the 
hadronic final state from moving fast in the B rest frame, and so the light-cone expansion 
which gives rise to the shape function is not relevant in this region of phase space.0 This 
is clear from the kinematics: the difference between the partonic and hadronic values of 
maximum is m%—ml ~ 2A nib, and nonperturbative effects are only important in a region 
of comparable size. For example, the most singular term in the OPE at order (AqcD/'^f))^ 
is of order (AqcD/'^ib)^ 5(1 — /ml). This contribution to the decay rate is not suppressed 
compared to the lowest order term in the OPE only if the spectrum is integrated over a 
small region of width Ag^ ~ Aqcd"^6 near the endpoint. This is the resonance region 
where only hadronic final states with masses mx ~ Aqcd can contribute, and the OPE 
is not expected to work anyway. In contrast, nonperturbative effects are important in the 

and m\ spectra in a parametrically much larger region where final states with masses 
m\ ~ Aqcd ^6 contribute.^ The better behavior of the spectrum than the Ei and m\ 
spectra is also reflected in the perturbation series. There are Sudakov double logarithms 
near the phase space boundaries in the Ei and m\ spectra, whereas there are only single 
logarithms in the q^ spectrum. 

The effect of smearing the q^ spectrum with the model distribution function in Eq. (1) is 



^The fact that the h quark distribution function is not relevant for large was pointed out in |22] 



in the context oi B ^ Xsi'^i decay and in |18| for semileptonic B decay. 

^Similar arguments also show that the light-cone distribution function is not relevant for small 
hadron energy. Ex, but it does enter for Ex near mb/2. If mb/2 — mc ^ Aqcd, then the constraint 
Ex < m£) in the B rest frame would also give a model independent determination of \Vub\- 
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FIG. 2. The lepton invariant mass spectrum to 0{as)- The meaning of the curves is the same 
as in Fig. |l|. 



Decay 


Width of region without 


Nonperturbative 


Fraction of 6 ^ n 


distribution 


charm background 


region near endpoint 


events included 




I^Ei = m%/2mB 


AEi ~ Aqcd 


~ 10% 




Arux = rnjy 


Arux ~ Aqcd mb 


80% 


dr/dg2 


Ag^ = 2mBmD — 


Aq^ ~ Aqcd "T-b 


~ 20% 



TABLE I. Comparison between the lepton energy, hadronic invariant mass, and lepton invari- 
ant mass spectra for the determination of |V^fc|. The region dominated by nonperturbative effects 
is parametrically smaller than the region without charm background only for the spectrum in 
the last row (viewing mjj ~ ~ Aqcd f^b)- The last column gives rough numbers corresponding 
to the plots in Figs. 1 and 2. 



illustrated in Fig. 0. In accord with our previous arguments, it is easily seen to be subleading 
over the region of interest. Table I compares qualitatively the utility of the lepton energy, 
the hadronic invariant mass, and the lepton invariant mass spectra for the determination 

of \Vub\- 

We now proceed to calculate the B — ^ Xud-v decay rate with lepton invariant mass above 
a given cutoff, working to a fixed order in the OPE (i.e., ignoring the light-cone distribution 
function which is irrelevant for our analysis). The lepton invariant mass spectrum including 
the leading perturbative and nonperturbative corrections is given by 

^ -(^ ' 2 (1 _ g2)2 (1 + 2g2) + 4 (3 - 45g^ + 30g6) 



Fq dg^ V 2m^y ^ / v ' ^ ' ' m 

+ ^^Xif) + ('^yp,Yif ) + ..., (5) 

where = q'^/ml, Pq = 11 — 2nf/3, and 



192 7r3 



(6) 



is the tree level b ^ u decay rate. The ellipses in Eq. (||) denote terms of order {Aqqb / 
and order terms not enhanced by (3q. The function Xlq^) is given analytically in Ref. p3| . 
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whereas Y{(f) was computed numerically in Ref. ||2^. The order l/m\ nonperturbative 



corrections were computed in Ref. |15|. The matrix element of the kinetic energy operator, 
Ai, only enters the (f spectrum in a very simple form, because the unit operator and the 
kinetic energy operator are related by reparameterization invariance [^]. Any quantity 
which can be written independent of the heavy quark velocity v must depend only on the 
combination (1 + Ai/2m^). The spectrum (and the total rate written in terms of m;,) 
are invariant under a redefinition of f , but, for example, the lepton energy spectrum is not 
since Ei = v ■ pe- ( Equivalent ly, the Ai term is a time-dilation effect, and hence is universal 
in any quantity that is independent of the rest frame of the B meson 

We shall compute the fraction of 5 — > X„£z/ events with > q^, F^Qq), as the relation 
between the total B — » Xuii^ decay rate and \Vub\ has been extensively discussed in the 
literature I^G],^, and is known including the full corrections |2^. After integrating the 



spectrum in Eq. (^, we can eliminate the dependence on the b quark mass in favor of the 



spin averaged meson mass ttib = (ruB + 3mB*)/4 — 5.313 GeV, following p9|. We find 



F{ql) = l-2Ql + 2Ql~Ql-^ (qI - 'iQl + 2g«) - ^ (g^ _ 7Q6 ^ ggs 

+ (Qo - 3Q^ + 2Qt) -^-^{Ql- ^Ql + Ql) 

+ ^xwS) + (^^y/3„yWS + ..., (7) 

where Qq = qo/rfiB- The functions X{Qq) and Y{Qq) can be calculated from X(g^) and 
y(g^). Converting to the physical B meson mass has introduced a strong dependence on 
the parameter A, the mass of the light degrees of freedom in the B meson. For = 
{niB- niof ^ 11.6 GeV^ we find F(11.6 GeV^) = 0.287 + 0.027a,(m;,) - 0.016a2(mb)/3o - 
0.20A/(1 GeV) - 0.02AV(1 GeV^) + 0.02Ai/(l GeV^J - 0.13A2/(1 GeV^) + . . . . The order 
a^A term is negligible and has been omitted. Using A = 0.4 GeV, Ai = —0.2 GeV^ [0 and 
ois{mb) = 0.22, we obtain F(11.6GeV^) = 0.186. There are several sources of uncertainties 
in the value for F. The perturbative uncertainties are negligible, as can be seen from the size 
of the Oi^alPo) contributions. At the present time there is a sizable uncertainty since A is 
not known accurately. In the future, a ±50 MeV error in A will result in a ±5% uncertainty 
in F. Finally, uncertainties from l/m\ operators can be estimated by varying the matrix 
elements of the dimension six operators within the range expected by dimensional analysis, 
as discussed in detail in [p!4| -[T6[|. This results in an additional ±4% uncertainty in F. We 



note that this is a somewhat ad hoc procedure, since there is no real way to quantify the 
theoretical error due to unknown higher order terms. Therefore, these estimates should be 
treated as nothing more than (hopefully) educated guesses. They do allow, however, for a 
consistent comparison of the uncertainties in different quantities. 

If q^ has to be chosen larger, then the uncertainties increase. For example, for q^ = 
15 GeV^ we obtain F(15 GeV^) = 0.158 + 0.024a,(m;,) - Q.Q12al{mb)(3o - 0.18A/(1 GeV) + 
0.0lAV(l_GeV2) + 0.02Ai/(l GeV^) - 0.13A2/(1 GeV^) + . . . ~ 0.067, using the previous 
values of A and Ai. The perturbative uncertainties are still negligible, while the uncertainty 
due to a ±50 MeV error in A and unknown dimension six matrix elements increase to ±14% 
and ±13%, respectively. (This uncertainty may be reduced using data on the rare decay 
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FIG. 3. (a) The fraction oi B Xud-v events with > Qq 



F{q^), to order e^LM and A^co/^^fe 
,3 



in the upsilon expansion. The shaded region is the uncertainty due to Aq^j-j/m^ terms, as discussed 
in the text. The dashed hne indicates the lower cut Qq = (m^ — m/j)^ ~ 11.6 GeV^, which 
corresponds to F = 0.178 it 0.012. (b) The estimated uncertainty in F^Qq) due to AQ(-;j-,/m^ terms 
as a percentage of F(gQ). 



B Xsi~^£~ in the large region, as discussed below.) 

Another possible method to compute F(go) uses the upsilon expansion p6|. By ex- 
pressing in terms of the T mass instead of rfiB, the dependence of -F(go) ^ ^''^^ 
is eliminated. Instead, the result is sensitive to unknown nonperturbative contributions to 
m-c- The uncertainty related to these effects can be systematically taken into account and 



has been estimated to be small One finds, 

, / B{B ^ XJiy)\^2^^2 1.6 ps\^^^ 
|VJ^(3.04±0.06±0.08)xlO-(^^^5^^^^-JLj . («) 

The errors explicitly shown in Eq. are the estimates of the perturbative and nonpertur- 
bative uncertainties in the upsilon expansion, respectively. 

For g2 = - mof we find ^(ll.eGeV^) = 0.168 + 0.016e + O.OUelj^^ - O.I7A2 + 
. . . ~ 0.178, where e = 1 shows the order in the upsilon expansion. This result is in 
good agreement with 0.186 obtained from Eq. (|^). The uncertainty due to A is absent in 
the upsilon expansion, however the size of the perturbative corrections has increased. The 
uncertainties due to 1/mf operators is estimated to be ±7%. For = 15GeV^, we obtain 
F(15GeV^) = 0.060 + O.Olle + O.Olle^LM - O-MAs + ... ^ 0.064, which is again in good 
agreement with 0.067 obtained earlier. For this value of Qq, the 1/ml uncertainties increase 
to ±21%. F{ql) calculated in the upsilon expansion is plotted in Fig. ^, where the shaded 
region shows our estimate of the uncertainty due to the 1/mf corrections. 

Concerning experimental considerations, measuring the q"^ spectrum requires reconstruc- 
tion of the neutrino four-momentum, just like measuring the hadronic invariant mass spec- 
trum. Imposing a lepton energy cut, which may be required for this technique, is not a 
problem. The constraint > (m^ — m^))^ automatically implies Ei > (m^ — rrioY /2mB ~ 
1.1 GeV in the B rest frame. Even if the Ei cut has to be slightly larger than this, the 
utility of our method will not be affected, but a dedicated calculation including the affects 
of arbitrary E^ and q^ cuts may be warranted. 
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If experimental resolution on the reconstruction of the neutrino momentum necessitates a 
significantly larger cut than = {mB — TnDf', then the uncertainties in the OPE calculation 
of Flq^) increase. In this case, it may instead be possible to obtain useful model independent 
information on the spectrum in the region > rn^^^s) — -'-"^•S GeV^ from the spectrum 
in the rare decay B Xsi^i~ , which may be measured in the upcoming Tevatron Run-II. 
There are four contributions to this decay rate, proportional to the combination of Wilson 
coefficients Cg, Cf^, CjCg, and C7. Cg is a g^-dependent effective coefficient which takes 
into account the contribution of the four-quark operators. Its g^-dependence yields negligible 
uncertainties if we use a mean Cg obtained by averaging it in the region 0.5 < g^ < 1 weighted 
with the b quark decay rate (1 — g^)^ (1 + 2g^). The resulting numerical values of the Wilson 
coefficients are Cg = 4.47 + 0.44 i, Ciq = —4.62, and C7 = —0.31, corresponding to the scale 
/i = rrih. In the g^ > rn^f^^s) region the Cy contribution is negligible, and the C^Cg term 
makes about a 20% contribution to the rate. For the Cg+Cf^ contributions nonperturbative 
effects are identical to those which occur in B ^ Xud-v decay, up to corrections suppressed 
by IC9 + C10I/IC9 — Ciol ~ 0.02. Therefore, the relation 

dr(5^X,£+£-)/dg2 \VM'' o? |C9|2 + |CioP + 12Re(C7C9)/(l + 2g2)' ^ ^ 

is expected to hold to a very good accuracy. There are several sources of corrections to 
this formula which need to be estimated: i) nonperturbative effects that enter the C-jCg 
term differently, ii) mass effects from the strange quark and muon, iii) higher cc resonance 
contributions m. B ^ Xgi'^i' , and iv) scale dependence. Of these, i) and ii) are expected to 
be small unless g^ is very close to m^. The effects of iii) have also been estimated to be at 



the few percent level [^, although these uncertainties are very hard to quantify and could 
be com 
a large 
by Eq. 



be comparable to the ±8% scale dependence |3T|] of the B Xsi~^i rate. Integrating over 
a large enough range of g^, q^ < q^ < Tn?B with m^^js) < Q'o ~ 17GeV^, the result implied 



BiB ^ XJ+i')\g2^^2 \VtsVt,\^ a2 |C9|2+|CioP + 12Re(C7C9)5(g2) 



(10) 



is expected to hold at the ~ 15% level. Here Qo = qo/rfiB, and -B(go) = 2/[3(l + Ql)] — 
4{K/fhB) (5o/[3(l + QlY] + . . . . For gg significantly above (m^ — ttidY, this formula may 
lead to a determination of \Vub\ with smaller theoretical uncertainty than the one using the 
OPE calculation of F(g^). 

In conclusion, we have shown that the g^ spectrum in inclusive semileptonic B — »• 
decay gives a model independent determination of \Vub\ with small theoretical uncertainty. 
Nonperturbative effects are only important in the resonance region, and play a parametri- 
cally suppressed role when dP/dg^ is integrated over g^ > (m^ — tjib,)^, which is required 
to eliminate the charm background. This is a qualitatively better situation than the extrac- 
tion of \Vub\ from the endpoint region of the lepton energy spectrum, or from the hadronic 
invariant mass spectrum. 
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